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A continuum theory of thermoelectric bodies and effective 
properties of thermoelectric composites 



Abstract 

We develop a continuum theory for thermoelectric bodies following the framework of continuum me- 
chanics and conforming to general principles of thermodynamics. For steady states, the governing 
equations for local fields arc intrinsically nonlinear. However, under conditions of small variations of 
electrochemical potential, temperature and their gradients, the governing equations can be reduced 
to a linear elliptic system and conveniently solved to determine local fields in thermoelectric bodies. 
The linear theory is further applied to predict effective properties of thermoelectric composites. In 
particular, explicit formula of effective properties arc obtained for simple microstructures of lami- 
nates and periodic E-inclusions, which imply useful design principles for engineering thermoelectric 
composites. 

1 Introduction 

Thermoelectric (TE) materials directly convert heat into electric energy. They are often used for 
power generation and refrigeration (Rowe 1999; Mahan 2001). Being all solid state and without 
moving parts, thermoelectric devices have the unique advantage of portability, silence, scalability, 
durability and reliability in extreme environment (DiSalvo 1999). Also, they arc readily customized 
and embedded into a large system to provide power or refrigeration. For example, thermoelectric 
materials have been used to power independent wireless systems (Leonov et al, 2008) and cool 
microelectronic chips (Boukai et al, 2008; Rowe 2006). Additional applications of thermoelectric 
devices include temperature and thermal energy sensors for precise temperature and thermal flux 
control (Riffat and Ma 2003), hybrid photovoltaic-thermoelectric energy systems which have an 
improved solar power conversion efficiency (Kraemer et al, 2007), etc. However, thermoelectric 
devices suffer from a critical disadvantage of low efficiency. 

Recent renaissance of research in thcrmoclcctrics has been focusing on improving thermoelectric 
efficiency, measured by a dimensionless figure of merit {ZTq). Equipped with modern technologies of 
nanofrabrication, tremendous amount of efforts have been devoted to explore cost efficient methods 
for manufacturing TE materials with high figure of merit, see, e.g., Yang et al. (2000), Venkata- 
subramanian et al. (2001), Yamashita et al. (2005), Ohta et al. (2007), reviews of Mahan (1998), 
Nolas et al. (2006), Snyder and Toberer (2008), Lan et al. (2010) and references therein. An 
important method is to consider TE nanocomposites which has demonstrated potential of improv- 
ing properties of TE materials including figure of merit ZTq, power factor, structural integrity 
and adaptivity, etc. Therefore, the quest for desirable TE materials would benefit greatly from a 
rigorous homogenization theory which can predict the macroscopic properties of TE nanocompos- 
ites and furnish design strategies for desired effective properties. Such a homogenization theory 
above all demands a continuum theory for thermoelectric bodies. Moreover, on the device level 
the geometry and boundary condition on the constituent thermoelectric body can significantly in- 
fluence local temperature and electrochemical potential on the body and hence the behavior of 
the device. The disparity between actual geometry and boundary condition in a functional device 
and ideal geometry and environment under which TE properties are measured may explain why 
the efficiency of TE devices is usually smaller than the efficiency of the TE material (Snyder and 
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Toberer 2008). It also calls into caution how to interpret experimental data when the geometries 
and boundary conditions of samples are significantly different from the ideal ones (Yamashita et al. 
2005; Yamashita and Odahara 2007). To quantitatively account for the influence of geometry and 
environment, we again need a self-consistent continuum theory capable of predicting all relevant 
local fields in thermoelectric bodies. 

As a continuum theory, the presented model of thermoelectric bodies shall consist of three 
components: kinematics, conservation laws and constitutive law, and this model shall conform to 
general principles of thermodynamics. Kinematically we describe the local thermodynamic state 
in a thermoelectric body by electrochemical potential and temperature (fi,T). The conservation 
laws can be formulated for energy flux j^, heat flux q, electric flux (or currents) je or entropy 
flux js, among which only two are independent. The constitutive law is usually postulated as a 
linear relation between fluxes and driving-forces which are typically gradient flelds derived from 
state variables {^,T). Such a constitutive law gives rise to a coefficient matrix L which may be 
measured by experiments or computed from microscopic theories (Ashcroft and Mermin 1976). The 
conformity to thermodynamic laws calls into question the validity of a postulated equation divjs = 
for steady states, used by, e.g., Bergman and Levy (1991) since 7 = divj^ can be identified as the 
rate of entropy production per unit volume for steady states (Gurtin et al., 2010). As a transport 
and irreversible process, a nontrivial steady-state of a thermoelectric body shall generate entropy 
and dissipate energy, i.e., 7 > 0. Moreover, there exist certain reciprocal relations in the coefficient 
matrix L as implied by Onsager's principle which is sometimes referred to as the "Fourth law of 
thermodynamics" (Wendt 1974). However, reciprocal relations hold only for "proper, conjugate" 
pairs of fiuxes and driving-forces, which gives rise to many confusions in the literature as for the 
"most convenient" pairs of fluxes and driving-forces. What make a bad situation worse lies in the 
temperature dependence of the coefficient matrix L. At a large temperature scale, the coefficient 
matrix L should inevitably depend on temperature T, which would make the macroscopic theory 
nonlinear and untractable; on the other hand we may assume the coefficient matrix L or any of 
its algebraic combinations with T and ji remain constant on a homogeneous thermoelectric body 
if variations of T and jj, are small. 

The above issues motivate us to scrutinize the underlying conditions or hypotheses for material 
properties, conservation laws, constitutive law, and important derived material properties such as 
figure of merit and power factor. We remark that essential elements of the proposed continuum 
theory are well-known to some extent in earlier publications, e.g., Donicnicali (1954), Callen (1948; 
1960), Barnard (1972), de Groot and Mazur (1984), Haase (1990), Milton (2002). Here our fo- 
cus is to develop the theory systematically following the framework of continuum mechanics. In 
particular, we clarify the conditions for applying self-consistently the proposed continuum theory: 
small variations of (//, T) and small driving-forces. The former condition is needed for assuming 
some algebraic combinations of the coefficient matrix and (/i, T) may be regarded as constant in a 
homogeneous TE body; the latter is needed for the postulated linear constitutive law and Onsager's 
reciprocal relations. These conditions are usually satisfied in practical applications of thermoelec- 
tric materials. If they were violated, a more complete model including temperature-dependence of 
material properties and nonlinear dependence of fluxes on driving-forces would be necessary. 

An important feature of the proposed theory is that the boundary value problem for determining 
local fields is a linear system of partial differential equations whose existence, uniqueness and 
stability have been thoroughly investigated, see e.g. Evans (1998). This is convenient since a 
nonlinear set of partial differential equations are rarely solvable, meaning reliable predictions are 
difficult and rare. Moreover, the established homogenization methods can be used for this linear 
system and we obtain the unit cell problem for predicting effective properties of TE composites. 
For simple microstructures of laminates and periodic E-inclusions, the unit cell problem admits 
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closed-form solutions, providing useful guide for designing TE composites. 

The paper is organized as follows. We begin with experimental observations and clarify condi- 
tions for defining conventional thermoelectric properties in Section 2.1. In Section 2.2 we formulate 
the conservation laws for energy flux j„ and electric flux je. In Section 2.3 we identity pairs of 
"proper, conjugate" fluxes and driving-forces to apply Onsager's reciprocal relations, and explain 
the rationale of assuming the coefficient matrix between chosen pairs of fluxes and driving-forces 
are constant on a homogeneous body. An alternative viewpoint of Onsager's reciprocal relations 
and positive-definiteness of the coefficient matrix is provided in Section 2.4. In Section 2.5 we 
formulate boundary value problems for various applied boundary conditions which can be used to 
determine local fields in TE bodies, and solve the problems for an infinite plate and circular or 
spherical shell in Section 3.1 and 3.2, respectively. We proceed to define the effective properties of 
thermoelectric composites by the homogenization theory in Section 4.1. Explicit formula of effective 
TE properties are obtained for composites with microstructures of simple laminates in Section 4.2 
and of periodic E-inclusions in Section 4.3. The applications of these closed-form formula to design 
desirable composites are explored in Section 4.4. We finally summarize in Section 5 and provide a 
derivation of Onsager's reciprocal relations for TE materials in Appendix A. 



2 A continuum model for thermoelectric bodies 
2.1 Experimental observations 

We first describe properties of thermoelectric materials observed in experiments. Consider a ho- 
mogeneous body C iR", where n = 2 or 3 is the dimension of space. The thermodynamic state 
of the body is described by electrochemical potential fi and absolute temperature T. Further, we 
view as primitive concepts such as internal energy density u, entropy density a, electric flux je, 
entropy flux j^, entropy supply r] (by external sources), and the identities 

q = Tj„ r = Tr?, (2.1) 

where q is heat flux and r is external heat supply inside the body (owing to radiation). Note that 
these quantities are pointwisely defined inside the continuum body Cl, permitted by the concept 
of quasi-closed and local equilibrium of subsystems in statistical mechanics (Landau and Lifshitz 
1999, page 2). Suppose that there is no external heat supply inside the body, i.e., r = 0; the only 
way of heat exchange with the environment is through the boundary of O by conduction. For a 
time-independent boundary condition on dQ, the body eventually evolves into a steady state for 
which the fluxes or the driving- forces, i.e., gradients of electrochemical potential and temperature, 
can be measured. Experimentally, the electric conductivity tensor a G IR"'^^, thermal conductivity 
tensor k £ IR'^^'^, Seebeck coefficient matrix s £ IBJ'^'^ and Peltier coefficient matrix (3 £ IR'^^'^ 
are introduced by the following observations: 

1. At zero temperature gradient, i.e., g = — (VT)-^ = 0, the electric flux je is given by the Ohm's 
Law: 

je = ere, e = -(V/i)^ on Q. (2.2) 

Note that the driving-force field e includes the electric field and chemical contribution from 
nonuniform electron concentration. 

2. At zero electric flux, i.e., je = (open circuit), the heat flux q is given by the Fourier's Law: 

q = Tj, = Kg, g = -{VTf on n. (2.3) 



4 



3. At zero electric flux, i.e., je = (open circuit), the driving- force field e = — (V//)^ for an 
electric charge is given by 



e = -sg, g = -(vr)^ on n. (2.4) 

The above coupling between the field e and temperature gradient g is called the Seebeck 
efi^ects. 

4. At zero temperature gradient, i.e., g = — (VT)-^ = 0, the heat flux q is given by 

q = Tj, = /3je on n. (2.5) 
The above coupling between the heat flux q and and electric flux je is called the Peltier effect. 

From (2.2)-(2.5), we infer the following phenomenological relation between fluxes and driving- 
forces: 



ji = Lifi, ji 



Li : = 



Tcr Tas 

/3(T {k + /3<TS) 



(2.6) 



where Li e iR2nx2n 

g, e are understood as n x 1 column vectors. In general the 
coefficient matrix Li shall depend on temperature T and independent of electrochemical potential /i. 
Further, by Thomson or Kelvin relations, we have Li = hj or equivalently, 

<T^ = <T, (3 = Ts^ , = K. 

The above symmetry can be regarded as a consequence of Onsager's principle of microscopic re- 
versibility (Onsager 1931 and Appendix). Further, from the Second law of thermodynamics, Li is 
positive definite, and henceforth for any nonzero vectors vi, V2 G iR", 

Tvi • <TVi + 2rv2 • crsvi + V2 • (k + rs'^crs)v2 > 0. (2.7) 

From a microscopic viewpoint the phenomenological linear relation (2.6) can be justified as the 
leading-order approximation when driving-forces on the carriers (i.e., electrons or holes) are weak 
(Ashcroft and Mermin 1970, chapter 13). In present context the driving-forces have two origins: 
electrochemical potential gradient and temperature gradient. To compare them self-consistently, 
we shall assume (cf., equation (13.43) of Ashcroft and Mermin 1970) 

(2.8) 

where ks = 1-38 x lO^^^J/X is the Boltzmann constant, e = 1.60 x lO^^^C is the charge of 
an electron, and Fq is the force such that the linear relation (2.6) is no longer applicable (when 
carriers are subject to driving forces at the order of Fq). This force Fq is a material's property and 
may be estimated from a pure electric measurement: Fq ~ gEq, where Eq is the critical electric 
field strength |e| such that the Ohm's Law j = ere is no longer sufficient to describe the flux-field 
relation. 
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2.2 Conservation laws 



Denote by p the free charge density and u the internal energy density. The conservation of electric 
charges implies 

divje = —dtp on 0. (2.9) 

Further, the internal energy density can be written as 

u = UQ + a,T + pp, (2.10) 

where uq is a constant independent of T, p, and Cy is the volumetric specific heat. Also, we identify 
energy flux as 

iu = Tis+pie, (2.11) 

and by conservation of energy, obtain 

divj^i = • VT + Tdivjs + je • V/x + ^divje = -dtu on O. (2.12) 



2.3 A constitutive model for thermoelectric materials 

There is a lot of confusion as to appropriate dependence of material properties (i.e., the coefficient 
matrix Li in (2.6)) on state variables {p,T) and how to justify the symmetry = Li by On- 
sager's principle. Prom a simple dimension analysis, we find that the symmetry Li = requires 
the dimensions of the two inner products between pairs of flux and driving-force, i.e., je • e and 



js • g, should have the same dimension. Further, let A G 
transformation of fluxes and driving-forces 



.2nx2n 



(ji,fi)^(j2,f2) = (Aji,A-^fi) 



be an invertible matrix. A 



(2.13) 



implies 



J2 



2^2, 



L2 = ALiA^ 



(2.14) 



Therefore, there is no unique choice of fluxes and driving-forces for describing the transport process. 
Moreover, we can switch the roles played by the flux and driving-forces and, e.g., refer to fs := 
[je,g/T] as the "driving-forces" and js := [e/T,js] as the "fluxes". Then by (2.6) we have 



Je 



33 



e/r 

is 



J3 



L,f- 



3I3; 



K 



(2.15) 



where the off-diagonal blocks in L3 are skew-symmetric instead of symmetric. Therefore, the 
Onsager's reciprocal relations do not always imply the symmetry of the coefficient matrix; the 
symmetric property of the coefficient matrix depends on the choice of driving-forces and fluxes 
(Coleman and Truesdell 1960). To remedy this issue, we subsequently require driving-forces be 
even functions of microscopic velocities of carriers whereas fluxes be odd functions of microscopic 
velocities of carriers, and hence the coefficient matrix would be symmetric by the Onsager's recip- 
rocal relations (Casimier 1945). 

It is clear that the inner products of pairs of driving-forces and fluxes remain invariant for 
transformations in (2.13) and (2.15): 



Ji 



•fi 



J2 



J3 • 13- 



(2.16) 
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Indeed, the fundamental thermodynamics for irreversible processes implies that the rate of entropy 
production per unit volume, denoted by 7, cannot depend on the choice of fluxes and driving-forces 
used to describe the irreversible process, and the conjugate pair of fluxes and driving-forces (j,f), 
by definition, satisfy (see Appendix A) 



7 = j • f- 

From Gurtin et al (2010, page 188) we have 

= dta + divjs, 

where a is the entropy density. By the First law of thermodynamics we have 

(9iu(x, t) = Tdta{-K, t) + ndtp{-x, t) , 
and by conservation laws (2.9)-(2.12), 

Tdivjs = -dtu - js • VT - jg • V/x iidtp. 
Inserting (2.20) into (2.18) and by (2.19) we obtain 

7 = (-je-v^-j,-vr)/r. 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



Therefore, (je,js) and (— V/z/T, — VT/T) are a conjugate pair of fluxes and driving- forces; the 
symmetry Li = in (2.6) follows from Onsager's principle of microscopic reversibility and the 
fact that pL and T are even functions of microscopic velocities of carriers. In Appendix A we outline 
a procedure of deriving Onsager's reciprocal relations for thermoelectric materials. 

For steady states, it is more convenient to use electric flux and energy flux (je, ju) since they are 
divergence free by conservation laws (2.9)-(2.12). To find their conjugate driving- forces, by (2.11) 
we notice the relation (I is the n x n identity matrix) 



J : = 



A 



A 



I ■ 



and hence, by (2.13) and (2.14), the conjugate driving-forces and coefficient matrix between the 
fluxes j (i.e., j = Lf) and driving-forces f are given by 



A" 



T 


--{Vpf/T 




"[v(-^/r)]^- 




_-(VT)^/T 




. [v(i/r)]^ _ 



To- fiTcr + TVs 



(2.22) 



respectively. 

As is well-known, an arbitrary additive constant in electrochemical potential fi should have no 
physical consequence. This makes the appearance of /x in the coefficient matrix L(/x, T) somewhat 
strange. Nevertheless, we recall that cr, s, k, are independent of /x (see discussions in Section 2.4). 
Upon straightforward calculations it can be shown that if {—p/T, 1/T) satisfy the conservation laws 
(2.9)-(2.12): 



div[L(^,r)f] 



'-dtp 


f = 


"v(-/x/r)" 


-dtu 




. v(i/r) _ 



(2.23) 
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then (— (/i + c)/T, 1/T) satisfy that for any constant c G IR, 



div[L(/x + c, r)f'] 



-dtp 
-dtu - cdtp 



Vi-p/T) - cV{l/T) 

v(i/r) 



(2.24) 



The above equation confirms that an additive constant in electrochemical potential p indeed has 
no physical consequence. 

From experiments or microscopic statistical mechanics, the coefficient matrix between fluxes 
and driving-forces in general depends on temperature. This dependence makes a boundary value 
problem based on conservation laws (2.9)-(2.12) nonlinear and untractable. This difficulty may 
be addressed by introducing an additional hypothesis that the deviation of the actual {p, T) from 
their equilibrium values (/io,7o) is small. For small bodies, this hypothesis is contained in the 
hypothesis of weak driving-forces which is required by the very linear relation between fluxes and 
driving-forces. Moreover, the Onsager's reciprocal relations require the system does not deviate too 
much from the equilibrium state which overlaps with the hypothesis of small variations of (//, T). 
Therefore, in typical situations it is reasonable to assume that the coefficient matrix is constant, 
independent of positions inside the homogeneous body fi, and given by 

L(^(x), r(x)) = L(/xo, To) Vx G n, (2.25) 

where (/io,2o) is the equilibrium electrochemical potential and temperature of the body when 
the body is closed and isolated. In practice we can choose Tq as the average temperature of the 
body and fio = 0. 

We remark that equally sound but not equivalent is to assume the coefficient matrix between 
another pair of conjugate fluxes and driving-forces, e.g., the matrix Li in (2.6), is uniform inside the 
homogeneous body. Compared with (2.25), this new assumption would imply different solutions 
to a specific problem, but the difference is presumably small and nonessential because of the same 
premises: weak driving- forces and small variations of (/i, T). However, as one will see below, 
equation (2.25) implies an easy-to-solve linear system of equations whereas Li being constant 
implies a nonlinear system (Yang et al., 2012). 

The temperature dependence of electric and thermal conductivities and Seebeck coefficient 
matrix implied by (2.25) is 

111 

o-oc-, Koc^, SOC-. (2.26) 

Ideally speaking, there could be thermoelectric materials with the above temperature dependence 
precisely over some temperature range. Even for such materials, the //-dependence of the coefficient 
matrix in (2.22) would still make (2.23) nonlinear for steady states. From this viewpoint, a rigorous 
continuum theory for thermoelectric bodies are intrinsically nonlinear. 

Finally we notice that the first and second of (2.26) may be justified by a microscopic theory 
for metals at a temperature much higher than the Debye temperature, see e.g. equations (26.48) 
and (13.58) of Ashcroft and Mermin (1976), whereas the justification of the last of (2.26) appears 
to be elusive, see comments in page 258 of Ashcroft and Mermin (1976). We stress here that the 
reasoning for (2.25) lies in the conditions of weak driving- forces and small variations of /x, T instead 
of an underlying microscopic theory, and that equation (2.25) can be applied to materials whose 
properties violate (2.26), e.g., semiconductors, as long as the assumed conditions are satisfied. 

2.4 An alternative viewpoint of the constitutive model 

The Second law implies that the rate of entropy production per unit volume 7 = f • Lf > for 
an irreversible process, which implies the positive-definiteness of the coefficient matrix L in (2.22). 
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Further, if the driving- forces f = 0, i.e., the electrochemical potential and temperature are uniform 
on the body, the body is in the equilibrium state and hence 7 = 0, attaining its minimum value. 
These observations motivate us to postulate the constitutive law between fluxes and driving-forces in 
irreversible processes by specifying the functional dependence of the entropy product rate per unit 
volume 7 on state variables and driving-forces. Below we derive the phenomenological linearity 
between fluxes and driving-forces, Onsager's reciprocal relations and positive definiteness of the 
coefficient matrix as consequences of this postulation and weak driving- forces. 

In analogy with a reversible process, e.g., the elastic theory of deformable bodies, and with the 
entropy product rate per unit volume 7 playing the role of internal energy density, we postulate 
that 7 is a smooth function of thermodynamic variables and driving-forces. The driving-forces, 
by definition, characterize the degree of deviation of a local subsystem from its equilibrium state. 
We consider a generic situation where the thermodynamic state is described by m variables: u = 
{ui, ■ ■ ■ ,Um) whose values would be uniform on the body in the equilibrium state and are even 
functions of microscopic velocity of carriers. Then the gradients F = Vu G iR™^" would be the 
appropriate quantities characterizing local deviation from local equilibrium state, and the above 
postulation can be expressed as 

7 = 7(u,F). 

Under the condition of weak driving- forces, i.e., |F| << 1, we have the Taylor expansion 

7(u, F) = 7o(u) + B(u) • F + F • C(u)F + • • • , (2.27) 

where 



7o(u) = 7(u,0), B(u) 



97(u,F) 



dF 



C(u) 



ia^7(u,F) 



p^o 2 dFdF 



(2.28) 

F=0 



Since 7(u, F) > if F 7^ (irreversible process) and 7(u, 0) = if F = (equilibrium state), we 
conclude that 



7o(u) = 0, B(u) = 0, C(u) is symmetric and positive definite. (2.29) 

From the definition, the proper, conjugate fluxes of the driving-forces F are the quantities J 
such that 



J-F = 7(u,F) = F-C(u)F, (2.30) 

where we have neglected higher order terms in (2.27). A natural solution to the above equation is 
given by 

{J)p, = {C)pigjiF)gj, i.e., J = C(u)F. (2.31) 

Therefore, we conclude that the postulation that the entropy production rate per unit volume 7 
is a smooth function of the driving-forces, to the leading order, implies the linear and reciprocal 
relations between the conjugate pairs of fluxes and driving- forces. 

When specified to thermoelectric bodies, if we choose u = (^, T) as the state variables F = Vu 
the driving- forces, we anticipate 7(u, F) and hence the associated fourth order tensor C(u) is 
independent of /x, which by (2.21) implies the coefficient matrix Li in (2.6) must be independent of 
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II. If we choose u = {ui,U2) = (— /x/T, 1/T) as the state variables, Vu as the driving-forces, then 
by (2.14), (2.21), and (2.22) we identify 



F = , J = I , J = C(u)F, 

where the associated fourth-order tensor C(u) G ^2xnx2xn g^ygj^ 

(C)iiij = T{a)ij, {C)ii2j = (C)2jij = fiT{a)ij +T^(crs)jj, 

(C)2i2i = [Ai^Tcr + /xr2(s^cr + crs) + T2(fc + rs^crs)]^^. (2.32) 

Further, for steady states, the conservation law (2.23) can be conveniently rewritten as 

div[C(u)Vu] = on J^, (2.33) 

where the divergence is taken over row vectors. 

We remark that the above viewpoint also facilitates deriving the implications of material sym- 
metry on material properties. For example, for isotropic materials we have 7(u, F) = 7(u, FR) 
for any rigid rotation R, which, together with (2.27), implies cr, k and s are isotropic matrix. A 
detailed discussion in a classic context can be found in Gurtin et al. (2010, §50). 

2.5 Boundary value problems for thermoelectric bodies in steady states 

We now consider a heterogeneous thermoelectric body 17 whose material properties may in general 
depends on position. If the boundary conditions applied to the body are time independent, the 
body would eventually evolve into a steady state. Suppose that the driving-forces to the carriers and 
variations of (;U, T) are small in the overall body, and we assume (2.25) such that the thermoelectric 
tensor of this body is given by C(uo, x), where uq = (— /io/^o, I/Tq) and (/^o, Tq) are electrochemical 
potential and temperature of the body in the equilibrium state. We remark that (/io,?o) are 
uniform on Vl even if the body is heterogeneous. Therefore, upon choosing /io = and setting 
u = (ui, U2) = (— ///r, 1/r), by (2.23) and (2.33) we infer that in steady states, the electrochemical 
potential and temperature fields satisfy 



Jdiv[Co(x)Vu] =0 on 0, 

[time-independent boundary conditions on dO., 



(2.34) 



where 



[Co(x)]iaj = ro[cr(x)]ij, [Co(x)]ii2j = [Co(x)]2jii = T^[a{-x.)s{-x)]ij, 
[Co(x)]2,2i = T2[fc(x) + ToS^cr(x)s(x)],„ (2.35) 

and the position-dependent cr(x),s(x), and k(x) reflect that the body may be heterogeneous. 
Further, let F C 9il be a sub-boundary of the body and 5, h be two given functions on F. 
Typical boundary conditions on F in applications include the following. 

1. Dirichlet boundary condition: 

T = g, fj, = h on F, (2.36) 
where g, h represents the temperature and electrochemical potential on F, respectively. 
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2. Neumann boundary condition: 



q-n = sr, je • n = /i on T, (2.37) 

where n is the outward normal on F, and g, h represents the outward heat flux and electric 
flux on r, respectively. 

3. Mixed boundary condition: 

Li(q-n,T)=5, L2{je ■ n, fi) = h on T, (2.38) 
where Li,L2 represent linear functions of two variables. 

Other types of boundary conditions are also allowed from the mathematical viewpoint but physically 
less common. 

A few remarks are in order here regarding transient states and sharp interfaces. Mahan (2000) 
suggested that for transient states the right hand sides of the conservation laws (2.9) and (2.12) 
should be given by 

dT 

dtp = 0, dtu = Cy — , 

where is the heat capacity per unit volume. The above equations, though commonly assumed 
when there is no coupling between heat flux and electric flux, are questionable for heterogeneous 
thermoelectric bodies since there could be accumulations of free charges in transient states and 
not all energy flow into a region would be converted into heat for a thermoelectric body. To close 
the system of equations and in analogy with electrostatics, we introduce the fleld of local electric 
displacement d and postulate an additional constitutive law: 

divd = p, d = — e(x)V/i, 

where e may be interpreted as the "permittivity" of the thermoelectric body. Then by (2.10) and 
conservation laws (2.9), (2.12) we obtain the following equations for transient states: 



div[Co(x)Vu(x,t)] 



9t[dive(x)V/x] 
-C,dtT + dt{pdiYe{^)Vfi) 



which appears to be a nonlinear system and difficult to solve. Moreover, it was suggested that 
sharp interfaces might play an important role in overall thermoelectric behaviors (Yamashita and 
Odahara 2007; Ohta et al 2007), which may be modeled by postulating a similar constitutive 
relation between fluxes and interfacial driving forces. Below we will not consider such interfacial 
effects for thermoelectric composites. 



3 Power factor and figure of merit of thermoelectric structures 

In this section we solve boundary value problems for simple structures of an inflnite plate and a 
spherical shell and extract some combinations of basic material properties, namely, power factor 
Pf and flgure of merit ZTq, for evaluating the capacity (per unit volume) and efficiency of power 
generation and refrigeration. 

For more complex structures or boundary conditions, we anticipate that the capacity and effi- 
ciency of power generation and refrigeration depend not only on material properties but also the 
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Figure 1: (a) An infinite thermoelectric plate of thickness L, electrochemical potential and temper- 
ature (Ti, — /Uo) on the bottom face and (T2,/xo) on the top face; (b) normalized maximum efficient 
^*/?ycarnot versus figure of merit ZTq. 



overall geometry and boundary condition, which can be predicted by solving the proposed bound- 
ary value problem described in Section 2.5. An interesting problem to investigate is whether there 
exist geometries and boundary conditions such that the overall capacity (per unit volume) and 
efficiency of power generation and refrigeration exceed those of an infinite plate described below, 
and if not, why. 

3.1 A thermoelectric plate 

Consider a homogeneous and generally anisotropic infinite thermoelectric plate of thickness L as 
shown in Fig. 1 (a). The temperature and electrochemical potential on top and bottom faces are 
maintained as 

T(z = 0) = ri, T{z = L)=T2, = 0) = -/io, fi{z = L) = fio- (3.1) 

A priori we shall assume the temperature difference 5T = T2 — Ti is small compared with the 
equilibrium temperature Tq, and hence 

1 1 2 1 1 (5r 

Also, we assume the variation of electrochemical potential is small and enforce the assumption (2.25). 
Let S iR" be the unit normal to the plate as shown in Fig. 1 (a) and denote by 

az = (C)mj(e^)i(e,)j/ro, a, = (C)ii2i(e,)i(e,)j/T2, 4 = {C)2i2j {eMez)j/Tl (3.3) 

which, by (2.35), can be identified as 

az=ez CTez, = • crse^, 4 = • (ft ros^<Ts)e^. (3.4) 

Let 

Ef = -^n, Pf = ^ (3.5) 
ct.k; az 
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be the efficiency factor and the power factor of the thermoelectric material. From (2.7), choosing 
vi = ae^ and V2 = be^ for any a,b £ IR, we see that 

az > 0, > 0, Ef < 1. 

For isotropic materials, the above parameters Ef and Pf are independent of directional vector 
and given by 



TnsV ZT( 



Ef = " ^ = = sV, (3.6) 

where ZTq = Tqs^it/k is the figure of merit of thermoelectric materials. By the first of (3.6) we 
have 

ZTo = (3.7) 
1 - Ef 

which remains valid for general anisotropic materials. Below we explore the physical meanings of 
ZTq and Pf. 

By the invariance of the problem for any translations within the plate, we infer that the solution 
to (2.34) and (3.1) satisfies 

Je = -To[— - —]ae. + ^0 - 

and hence the normal electric and energy fluxes along e^-direction is given by 

Tnllna, ,1 1 ^ T^oiz ,1 1 X "^IJ-of^z STUz 

Ju-3u-ez- ^ + ^ [^^ ^ ^ . 

Therefore, the rate of work per unit volume done by the plate, the heat flux into the plate from 
the bottom and the heat flux out of the plate from the top is given by 

W = 2jel2o/L, qi = ju + Je/Ho, 92 = ju " je/^O, (3-9) 

respectively. Let 

5Taz 



be a dimensionless number for future convenience. By (3.8), (3.5) and (3.10) we find that 



(3.10) 



Je/xo~ —{l + p), ju^ ir~6f^'^ 
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Power generation 

Without loss of generality we assume /j-q > 0. If qi > q2, i.e., je > p < —1, the thermoelectric 
plate converts heat into electric energy. The electric power generated per unit volume is given by 

2 2 

W := 2fioje/L = -j^{2^iQaz + no5Taz). 

Therefore, for fixed applied temperature difference 5T the maximum power generated per unit 
volume is given by 

We remark that power factor Pj is more important than efficiency (figure of merit ZTq) for a power 
generation system where heat sources are abundant and cost is negligible (Rowe 1999; Narducci 
2011). 

Further, the efficiency of power generation is given by 

qi - q2 2fioje 2 

^ max{|gi|,|g2|} \ju\ + fioje + + p)| + l ' 

Upon maximizing the efficiency against p < — 1, we find 

To ^1 + ZTo + 1 + <5r(Vl + ZTo - l)/(2ro) ' ^' ^ 

which is achieved at p = —1 — + ZTq. We remark that since 6T/Tq « 1, the the above 

formula is well approximated by 

,_6T ^l + Zn-l 

which agrees with the optimal efficiency predicted by loffe (1957, page 40). 
Refr igerat ion 

Without loss of generality we assume (5T>0. Ifgi>0, i.e., 

2 ST 
qi = -jWzpI + in + 6T/2)a,po + ^4] > 0, (3.14) 

the thermoelectric plate extracts heat from the low temperature side and ejects heat into the high 
temperature side. Since o"^ > 0, the above inequality is possible if and only if 

{To + 6T/2fal-2a.K'JT>0 (1 + |^)2 _ ^ ^ > q_ (3_^5) 

Therefore, 

0<6T< := 2ro ^l±f^~| , (3.16) 

V 1 + ^-^0 + 1 

which is consistent with the prediction of Goldsmid (2010, page 11) when ZTq is small and such 
that 6T/Tq « 1. We also remark that other possibilities implied by (3.15) are omitted for 
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Figure 2: (a) A thermoelectric shell with prescribed temperatures and electrochemical potentials 
on the two faces; (b) the geometric factor Kn{R2/ Ri) versus R2/R1 in (3.23). 

contradicting our assumption 6T/Tq << 1. For the same reason, STj^ax, the maximum temperature 
difference possibly maintained by a thermoelectric cooler, cannot serve as a good measurement of 
ZTo if 6T max/To is not a small number. In this case, a more detailed analysis accounting for the 
temperature dependence of material properties is necessary to predict (JTmax- 
Further, by (3.11) the cooling efficiency is given by 

r]c- 



-2j,fio 2(1 + p) 

where p > —1 since jeP^o < 0, meaning work is done to the thermoelectric plate. Upon maximizing 
rjc against p > —1 we find 



ric - gj, 



yi + ZTq - 1 _ 6T 
VI + ZT0 + I 2To 



at p = -l + 1/Vl + ZTo, (3.17) 



which agrees with the prediction of Goldsmid (2010, page 10) when ST/Tq « 1. 

Recall that the ideal efficiency for a reversible system working between two heat reservoirs at 
temperature Ti and T2 is given by the Carnot efficiency: 

f?carnot ~ TfT for power generation f« — ^ for refrigeration, 
To oT 

where the temperature difference ST = \T2 — Ti\ << Tq = -^{Ti +T2). Therefore, a common factor 
between the optimal efficiency and Carnot efficiency can be extracted out from (3.13) and (3.17): 

7]* VI + ZT0-I 



r/carnot VT+~ZTo + 1 ' 

whose functional dependence on the figure of merit ZTq are shown in Fig. 1 (b). 
3.2 A thermoelectric shell 

In applications it can be advantageous to use geometries other than plates/films for structural 
adaptivity, ease of thermal insulation or conduction, etc. In this section we explore the geometric 
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effect by considering spherical shells in two dimensions and three dimensions as shown in Fig. 2 (a). 
For simplicity, we assume the thermoelectric material is isotropic with electric conductivity o", 
Seekback coefficient s and thermal conductivity k, and the boundary conditions are given by 



T{r = Ri)=Ti, T{r = R2)=T2, /x(r = = -//q, n{r = R2) = ^q. 



(3.18) 



Again we assume the temperature difference 6T = T2 — Ti is small compared with the equilibrium 
temperature Tq, and hence (3.2) holds. Below we derive the power factor and efficiency of this 
thermoelectric shell. 

Since the material is isotropic, equation (2.34) can be written as 

{Aui + TqsAu2 = on 17, 

saAui + (k + Tocrs^)An2 = on 0. 

From the symmetry we seek a solution to the above equations which can be written as u = 
{ui{r),U2{r)). Direct calculations show that 



(ui + Tosua) = 0, yreiRi,R2), 

{saui + {K + Tos^a)u2) = 0, Vr G {Ri,R2). 



1 


^ n- 




n— 


r 

^ dr 


dr 


1 


d n- 


i± 


n— 


r 

^ dr 


dr 



The general solution to the above equations is given as follows: 



if n 



ui = -Ax log(r) + Ao, 
U2 = --Bilog(r) + Bo; 



if n = 3, 



T 



+ A 



0, 



U2 = ^+Bq, 

where Ai,Aq,Bi and Bq are constants. By boundary conditions (3.18) and (3.2) we find 

^1 - Tolog(i?2/Ri)' ^l'- 

^0- rir2iog(i?2/iJi) ' o ° 

if n = 3, < 

Bi « ^,1 ,^ Bi 

Bo 



if n 



Bn 



, _ ST 

'1 ^ r2iog(ij2/Hi)' 



ro(iJ2-Ri)' 

(RiT2+Tifi2)A'0 
(R2-Ri)T^T2 

R2R1ST 
iR2-Ri)Ti' 



(3.19) 



. >0 - TiT2log(R2/Ri) 

Therefore, the electric flux and energy flux in radian direction are given by 



T1R2-R1T2 

TiT2{R2-Ri)- 



ije{Ri) = -{Ai + TosB,)ToaRl-'\ iju{Ri) = - 
\je{R2) = -{Ai + TosB,)ToaRl-^, \ju{R2) = - 

and hence the heat flux in radian direction is given by 



[Aias + Bi{k + Toas^)]T^R\-'', 
[Aias + Bi{k + Toas^)]T^Rl-'', 



(3.20) 



q{Ri) = ju{Rl) + ^oJe(-Ri), 

^(-^2) = i«(^2) - ^^0je{R2)■ 



(3.21) 
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Power generation 

Without loss of generality we assume fiQ > 0. If je > 0, i.e., 

Ai + TqsBi < 2fio + s6T < 0, (3.22) 

the thermoelectric shell converts heat into electric energy. The electric power generated per unit 
volume is given by 

where Kn{R2/Ri) is a dimensionless geometric factor given by 

Kn{x) 



if ^ = 2- 

(a::+l)logx ' 

if n = 3. 



3x 



Therefore, for fixed applied temperature difference 6T the maximum power generated per unit 
volume is given by 

= —^R^^^ ^0 = -^' (3.23) 

where Pj = cjs^ is the power factor of the thermoelectric material. Compared with a plate, the 
additional factor X„(|^) in the above equation is plotted in Fig. 2 (b), from which we observe that 
the influence is negligible for thin shells, i.e., R2 — Ri « Ri- 
Further, the efficiency of the conversion is given by 

^ 2/ioje(i^i) ^ 2(j(-2/ig - fiQsdT) 

\ju{Ri)\ + We(-Ri) |2as/xoro + {k + Toas^)6T\ - cj(2/x2 + ^os6T) ' 

Upon maximizing the above efficiency again satisfying (3.22), we find the maximum efficiency is 
again given by (3.12) with ZTq = T^s^a/ k and achieved by //q = —STs^Jl + ZTq{-\/1 + ZTq — l)/2. 

Refrigeration 

Without loss of generality we assume 5T > 0. Ifg(i?i) > 0, i.e., 

- 2asfioTo -{k + Toas^)6T - a{2fil + hqsST) > 0, (3.24) 

the thermoelectric plate extracts heat from the low temperature side and ejects heat into the high 
temperature side. Since a > 0, the above inequality is possible if and only if 

(To + (5r/2) V^s^ - 2a{K + Toas^)6T > 0. (3.25) 

Therefore, 



0<5T< 5r„,,. := 2To ^l±|^ | , (3.26) 

which is the same as (3.16). Further, by (3.20)-(3.21) the cooling efficiency is given by 

_ qjRi) 

and upon maximizing tJc against /^o satisfying (3.24) we find the maximum cooling efficiency is 
again given by (3.17). 
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4 Heterogeneous thermoelectric media 



(a) 




c c c 



II 



Figure 3: (a) A representive volume element of multiphase composites with the rth phase occupying 
^r] (b) a two-phase laminate composite of Ci and C2. 



4.1 Effective properties of thermoelectric composites 

We now consider a heterogeneous medium of N individual thermoelectric phases as shown in 
Fig. 3(a). Let Y = U^^ri^ be a representitive volume element or unit cell, and Q..,- C 1^ be the 
domain occupied by the rth phase (r = 1, • • • , N). Suppose weak driving- forces and small variations 
of (/i, T) in the overall composite body and enforce the assumption (2.25) for each individual phase. 
Note that the electrochemical potential ^uq is constant in the equilibrium state even if the body is 
heterogeneous. As for a homogeneous body, we choose /Uq = and set u = (— /i/T, 1/T). From the 
homogenization theory (Milton 2002), if the overall composite body is much larger than the unit 
cell Y , then the thermoelectric behavior of this composite body can be regarded as an "effectively" 
homogeneous body with an effective thermoelectric tensor C^. This effective tensor can be 
determined by solving the unit cell problem: 



(4.1) 



J div[C(x)(Vv + F)] = on y, 
1 V is periodic on dY, 

where v = u - Fx, F = Vu G m?^"^ is the average of the gradient of u in a unit cell, 

C(x) = C^ ifxea,, r = l, 

(Cr)ljlj = TQ{(^r)ij-, {Cr)li2j = {Cr)2jli = T^{arSr)ij, (4.2) 
(Cr)2j2j = [T'o(«^r + TqS^ (T rS,j.)]ij , 

Tq is the equilibrium temperature, and cr^, s^., k,. (r = 1, • • • , A^) is the electric conductivity tensor, 
Seeback coefficient matrix, and thermal conductivity tensor of the rth phase, respectively. 

By definition, the effective tensor C^, mapping local average Vu to local average fluxes, is such 
that 



j- ^ = j- C(x)Vu = C^y" Vu = C^F VFeiR^xn^ 



(4.3) 



where -^y for a domain V denotes the average of the integrand over V . From (4.3), we can 
alternately define the effective tensor by the quadratic form 



F • C"F = F • C (x) Vu VF e iR^""", 



(4.4) 
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which, in account of (4.1) and the divergence theorem, is equivalent to 



F-CT = niin|^ (Vv + F) • C (x) (Vv + F) | VF G iR^x'^. 

Here the admissible space W for v includes the completion of all periodic continuously differentiable 
functions from Y to IR^ in a suitable norm (i.e., W^''^). The above variational formulation can be 
regarded as reminiscent of the principle of minimum entropy production of Prigozhin (1962) and will 
be useful for deriving bounds on the effective tensor C^. We will not address bounds in this paper. 
An excellent account of bounds on thermoelectric composites or more generally, elliptic systems 
with unequal dimensions of domain and range, has been presented in Bergman and Levy (1991) 
and Bergman and Fel (1999), where they proved that the figure of merit of two-phase composites 
is bounded from above by the larger one of the constituent phases if the constituent phases and 
the composites are isotropic. 

Below we calculate the effective thermoelectric tensor of two-phase composites with microstruc- 
tures of simple laminates and periodic E-inclusions. These solutions are rigorous and closed-form, 
giving useful guide on designing thermoelectric composites for a variety of applications. 

4.2 Effective thermoelectric properties of simple laminates 

As illustrated in Fig. 3 (b), we now consider a two-phase laminate composite of material property 
tensor Cr, volume fraction 9^ {r = 1,2), and interfacial normal n. It is well known that the cell 
problem (4.1) can be solved explicitly for laminates, see e.g. Milton (2002, page 159 and references 
therein). Instead of quoting the final formula, we outline below the calculations for the reader's 
convenience. 

Let u = Fx -|- V, V be a solution to the cell problem (4.1). Clearly the gradient field Vu and 
the fluxes C(x)Vu are constant within each laminate, and are denoted by F^ = Vu on il.^ and 
Jr = C^Vu on Oj. (r = 1,2), respectively. The potentials u being continuous across interfaces 
implies that Fi — F2 = a (8) n for some vector a G IR^. Also, we have ^iFi -|- 62F2 = Vu = F, 
and henceforth 



Fi = F + e2a®n, F2 = F-0ia«)n. (4.5) 
Further, the conservation laws across interfaces imply 

(Ji-J2)n=(CiFi-C2F2)n = 0. (4.6) 
Inserting (4.5) into (4.6) we obtain 

[(Ci - C2)F]n + [(02C1 + 0iC2)a ® n]n = 0. (4.7) 
Moreover, by definition we have 

J = ^C(x)Vu = ^iCiFi + e2C2F2 = (^iCi + ^2C2)F + 6162(01 - C2)a ^ n =: C^F. (4.8) 

For future convenience, we introduce notations: 

Cg = 62C1 + 61C2, Cg = 61C1 + 6202- 
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Also, let N be the 2x2 inverse matrix of {Ce)piqj{'n)i{ii)i {p,q = 1,2), and S : iR^^" 
a symmetric tensor with its components given by 

(S)pigi = (N)pg(n)i(n)j. 

Upon solving (4.7) for the vector a and inserting it into (4.8) we find 

= Cg — 0i^2(Ci — C2)S(Ci — C2). 

If the constituent phases are isotropic in the sense that 

fTf> — CT^I^ S — S^Ij — /S^j'I ^T" — 1^ ^)) 



iR2xn be 



(4.9) 



or equivalently (cf. (4.2)), 

{Cr^piqj = (A7-)pq(5jj, Ay. 

then equation (4.9) implies the effective tensor has its components given by 



7q (TfSy 



Jptqj 



{Ag)pq6ij 



ei92{n)i{n)j{Ai - A2)pp'(Ag ^)pv(Ai - A; 



2)q'q, 



(4.10) 



(4.11) 



where Ag = 9iAi + ^2A2 and Ag = ^2Ai + ^iA2. For applied boundary conditions as shown in 
Fig. 1(a), we denote relevant effective material properties by 

Tq az Tq 

where az,az,Kz are given by (3.3) with the tensor C replaced by in (4.11). If n = e^, i.e., 
the lamination direction n is parallel to the applied temperature and electrochemical potential 
gradients, by (4.11) we find 



hA^^+92A^Y'- 



(4.12) 



If 62 • n = 0, i.e., the lamination direction n is perpendicular to the applied temperature and 
electrochemical potential gradients, by (4.11) we find 



7iAi + 02A2. 



(4.13) 



We remark that, though equations (4.12) and (4.13) resemble the classic mixture rules for 
simple laminates, equation (4.12) (resp. (4.13)) implies that the effective electric (resp. thermal) 
conductivity of thermoelectric composites is in general not given by the familiar mixture rule: 
(^i/cTi + ^2/c2)~^ (resp. 9iKi + 92K2), for the coupling between electric and thermal fluxes or the 
nonzero off-diagonal components in Ai and A2. 

4.3 Closed-form solutions for periodic E-inclusions 

Particulate heterogenous media are frequently encountered for which the estimates (4.9) and (4.11) 
are not suitable for their obvious dissimilarity in microstructure. It is therefore desirable to rigor- 
ously solve the unit cell problem (4.1) for microstructures resembling a particulate composite and 
obtain its effective material properties. This is usually achieved by using the Eshelby's solution 
for ellipsoidal inclusions in the dilute limit where the interactions between various particles are 
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neglected (Eshelby 1957; Mura 1987). For finite volume fractions and to account for mutual inter- 
actions between particles, we consider recently found geometric shapes called periodic E-inclusions, 
for which the unit cell problem (4.1) can be solved in a closed- form as for ellipsoidal inclusions in 
the dilute limit. 

We consider two-phase thermoelectric composites with their thermoelectric tensors given by 

C(x) = Cr if X G = 1,2, 

where is referred to as the inclusion (occupied by a second phase thermoelectric particle), r^i 
represents the matrix (continuous) thermoelectric phase, and Cr {r = 1,2) are of form (4.10). We 
remark that many but not all anisotropic thermoelectric tensors can be transformed into this form 
under linear transformations. 

Below we solve explicitly the unit cell problem (4.1) for 0.2 being a periodic E-inclusion following 
the well-known Eshelby equivalent inclusion method. First, we recall that a periodic E-inclusions, 
by definition, is such that the overdetermined problem (Liu et al. 2007; Liu et al. preprint) 

v2(/> = - xn^ on y, 

VV0 = -(1 - 0)Q on O2, (4.14) 

periodic boundary conditions on dY, 

admits a solution, where = \Vt2\/\Y\ is the volume fraction of the inclusion, and the shape matrix 
Q is necessarily symmetric positive semi-definite with Tr(Q) = 1. Then, by Fourier analysis or 
Green's functions, we find solutions to the following homogeneous counterpart of (4.1) 



V • [CiVv + T,*xn.,] =0 on y, 

periodic boundary conditions on dY 



(4.15) 



satisfy that, for any G IR^^^, the gradient field Vv is uniform on Q2 and is given by 

Vv = -(l-0)(Ai)-iS*Q=: -(1-0)RS* on Jls, (4.16) 

where the components of the tensor R : IR^^^ — )• IR^^^^ are given by 

(R)pig, = (Ai);^^ (Q)i, {p,q = l,2;i,j = !,■■■ ,n). (4.17) 

We now consider the inhomogeneous problem (4.1). Based on the equivalent inclusion method we 
observe that the solution to (4.1) is identical to that of (4.15) if the average applied field F for 
(4.1) and the "eigenstress" for (4.15) are related by (Liu et al. preprint) 

ACF = (1 - e)ACRY;* -S* = [(1 - 9)aCR - n]5]*, (4.18) 

where AC = Ci - C2 and II : iR^^" IR^^'^ is the identity mapping. 

To calculate the effective tensor of the composite, by (4.4) and (4.16) we find that the effective 
tensor satisfies 

¥-C^V = j- F • (Ci - ACxqJ (Vu + F) dx 

= F • CiF - OY ■ aC[-(1 - 6i)R5]* + F]. 
By (4.18) we rewrite the above equation as 

F-C^F = F-CiF + 6lF-S;*. 
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Figure 4: (a) Contours of normalized effective power factors P^/Pp and (b) contours of normalized 
effective power factors Z^Tq. The horizontal and vertical axes represent {a, k) (cf., (4.21)), and s is 
such that the figure of merit of the a-phase is the same as the i?-phase. The composites consist of 
a-phase periodic E-inclusions in i?-phase matrix with isotropic shape matrix and volume fraction 
Oa = 0.5. 

Further, it can be shown that the tensor (1 — 0)aCR — II is invertible for generic cases and the 
above equation implies 

= Ci -^)ACR-n]-^AC, (4.19) 

which is our closed- form formula of effective tensors for two-phase thermoelectric composites. 

We remark that formula (4.19) is a rigorous prediction to the effective properties of thermo- 
electric composites with microstructures of periodic ^'-inclusions and there is no phenomenological 
parameters in (4.19). Further, the shape matrix Q can be any positive semi-definite matrix with 
Tr(Q) = 1 and the volume fraction 9 can be any number in (0, 1), which allows for modeling 
composites with a variety of textures and volume fractions of inhomogeneities. In particular, if 
the shape matrix Q is a rank-one matrix n n for a unit vector n, the inclusion degenerates to a 
laminate perpendicular to n and (4.19) recovers the formula (4.11) for simple laminated composites 
whereas if Q = I/n we have 

{C%,gj = {A%g6ij, A^ = Ai + 0[(1-0)aAA^V"-I]"'aA. (4.20) 
4.4 Applications to the design of thermoelectric composites 

The closed-form formula (4.19) are useful for engineering the efficiency and power factor of thermo- 
electric composites. To this end, we choose p-type Bismuth-telluride (Bio.25 Tio.75 Tes doped with 
6 wt% excess Te) as the base material (referred to as the i?-phase). The thermoelectric properties 
of this p-type Bismuth-telluride are listed below (Yamashita and Odahara 2007): 

ab = 0.326 X 10^(^^m)-^ Sb = 245.0/iF/i^, Kb = 0.559W/mK, 

by which the power factor and figure of merit can be easily calculated and are given by 

P^f = 2.0 X W-'^^W/mK'^, ZbTo = 1.04. 
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Suppose the base material be mixed with a second isotropic thermoelectric material (referred to as 
the A-phase) with 



0"a = crab, Sa = SSb, Ka = KKb, (4-21) 

where a, s, k are, respectively, dimensionless scaling factors of electric conductivity, Seebeck 
coefficient, and thermal conductivity. Below we assume a family of shape matrices Q parameterized 
by a scalar uj >\: 



Q(. 







(4.22) 



We remark that if a; — ?• and > 0, the periodic E- inclusion approaches to a fiber 
and if a; ^ oo and ^ > 0, the periodic E- inclusion approaches to a laminate in G_2-direction. From 
(4.19) we can write the effective power factor and figure of merit along e^-direction as functions of 
Cl), 6, (7, s, k: 

= Pf{u;,e,a,s,k), 
Z^To = Z^To{uj,e,a,s,k). 

Below we investigate numerically how the effective power factor and figure of merit depend on these 
parameters. 

Since it is hard to improve the figure of merit of thermoelectric materials, we first vary a, k 
between (10~^, 10'^) and choose s = ((t/k)^/^ such that the figure of merit of A-phase is the same 
as that of i?-phase. We assume the microstructure of composites are A-phase particles of the shape 
of periodic E-inclusion embedded in i?-phase matrix. The volume fraction of inclusions is fixed at 
0.5 and the shape matrix is given by (4.22) with u = 1. Figure 4 (a) shows contours of normalized 
effective power factors P^/P!^ whereas Fig. 4 (b) shows contours of effective figure of merit Z^Tq. 
From these figures we observe that for fixed thermal conductivity there exists an optimal electric 
conductivity for maximizing the effective power factor and vice versa. However, the power factor 
cannot be significantly improved (within five fold). Also, as shown in Bergman and Levy (1991) 
the effective figure of merit cannot exceed that of the constituent phases. Moreover, if electric and 
thermal conductivity are uniformly scaled, i.e., a = k, the effective figure of merit is the same as 
that of the constituent phases. This remains true even if s = k = 0, which implies that porous 
media of high figure of merit materials would likely keep the property of high figure of merit (Inoue 
et al 2008). Additionally, we remark that similar qualitative results are observed when volume 
fraction 6 and shape factor uj are varied and when inclusion material and matrix material are 
interchanged. 

Next we consider composites of copper (A-phase) periodic E-inclusions embedded in B-phase 
matrix. The scaling factors are given by o" = 2.0 x 10^, s = 7.8 x 10^^ and k = 7.17 x lO'^ for 
copper (Yamashita and Odahara 2007). Figure 5 (a) & (b) show respectively contours of effective 
power factor Pj/Pj and effective figure of merit Z'^Tq, where the horizontal and vertical axes 
represent volume fraction and shape factor {6,uj). Qualitatively, the shape of periodic E-inclusion 
is like a needle (disk) when oj « 1 {uj >> 1). From Fig. 5 (a) we observe that the effective 
power factor increases as the volume fraction of j4-phase increases. In contrast, there exists an 
optimal value of uj for maximizing power factor. As volume fraction 9 increases from 0.1 to 0.9, 
the optimal shape gradually changes from needle-like shapes to disk-like shapes. Also, the power 
factor can be significantly improved (up to 100 fold). From Fig. 5 (b) we observe that disk- like 
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Figure 5: (a) Contours of normalized effective power factors Pj /^/ ) ^''^'^ (b) contours of normalized 
effective power factors Z^Tq. The horizontal and vertical axes represent volume fraction and shape 
factor {9a, uj) (cf., (4.22)), and the composites consist of vl-phase periodic E-inclusions in i?-phase 
matrix. The material properties of A-phase takes values of copper. 

shapes in general have higher figure of merit than needle-like shapes. This dependence is however 
weak in the sense that the change of effective figure of merit is within 10 ~ 20% even if uj varies 
between (0.1, 10^). Therefore, the trade-off between power factor and figure of merit often makes 
shapes with significantly improved power factor but a sightly lower figure of merit more desirable. 
Figure 5 (a) & (b) provide valuable information for choosing optimal volume fractions and shapes 
of copper particles to fulfill practical requirements. Finally we remark that interchanging inclusion 
and matrix phase materials yields undesirable lower power factor and figure of merit. 

5 Summary and discussion 

We have developed systematically a continuum theory for thermoelectric bodies following the gen- 
eral framework of continuum mechanics and conforming to basic thermodynamic laws. Under 
reasonable assumptions of small variations of electrochemical potential and temperature, the re- 
sulting boundary value problem is a linear system of partial differential equations which can be 
conveniently solved to determine all relevant local fiends. We have also applied the theory to pre- 
dict the effective properties of thermoelectric composites. In particular, explicit formula of the 
effective thermoelectric tensors have been calculated for simple microstructures of laminates and 
periodic E-inclusions. These explicit formula facilitate investigations of how practically important 
material properties such as power factor and figure of merit depend on the volume fractions, shapes 
of inclusions and material properties of constituent phases. 

Appendix A: Onsager's reciprocal relations and their applications 
to thermoelectric materials 

The original Onsager's derivation of the reciprocal relations of thermoelectric materials is somewhat 
elusive to researchers in the community of continuum mechanics. Following Casimier (1945) we 
outline below the premises for deriving the reciprocal relations for thermoelectric materials whereas 
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detailed derivations can be found in Onsager (1931) and Casimier (1945). 



General systems in irreversible processes 

We consider a closed isolated system. Assume the system is described by a set of macroscopic 
thermodynamic variables a* {i = 1, • • • , A^) whose values, without loss of generality, are assumed 
to be zero in the equilibrium state. Although the system is in the equilibrium state, the system 
can fluctuate, which is characterized by time-dependent random variables a(i) = (a^, • • • , o'^)^. 
The derivation of Onsager's reciprocal relations requires the following premises: 

1. The entropy in a state a deviating from the equilibrium state, to the leading order, is given 

by 

5(a) = S-o - • Sa, (A.l) 

where the constant N x N symmetric matrix S is a property of the system, and positive 
definite since entropy is necessarily an maximum in the equilibrium state. From the definition 
of entropy (Kittel and Kroemer 1980, page 40), the number of microscopic admissible states 
of the system is given by g{a.) oc exp(— g^a • Sa), and the probability density function is 
given by 

T^(a) = iexp(-^a.Sa), 



where C = ^{27r)^k^/ det S is a normalization constant such that J^jv VF(a)(ia = 1. Let 
b = Sa. Direct integration shows that 

E[h(g>ai]= I bOaVF(a)o?a = feijl, (A.2) 

where I is the N x N identity matrix. 

2. Also, we assume the ergodicity of all admissible microscopic states such that the expectation 
of a quantity /(a) is exactly equivalent to a time averaging: 

i?[/(a(t))]= hm 1 r/(a(t))dt. (A.3) 

T->+oo I Jq 

3. Let 

a(t,T) = E[a(t + r)|a(t)] 

be the conditional expectation characterizing the probability of finding the system in states 
a{t + r) at time t + t while knowing the system is in state a(t) at time t. We consider only 
state variables a which are even functions of velocity (Casimir 1945). Then the microscopic 
reversibility implies 

a(t,T) = a(t,-r). (A.4) 

Therefore, 

E[a{t) (g> a(t + r)] = E[a{t) a{t, r)] = E[a{t) (g> a(t, -r)] = E[a{t + t) (g> a{t)], (A.5) 

where the first equality follows from the general probability theory (see e.g. Evans 2002, page 
31), the second equality follows from (A.4), and the last equality follows from (A.3). 
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4. If for irreversible nonequilibrium processes the evolution of the system follows the macroscop- 
ical kinetic law for a constant N x N matrix M: 

a(t) = Mb(t), (A.6) 

then the fluctuation (in the equilibrium state) follows the above law as well in the sense that 

E[si{t + r) - a(t)|a(t)] = TMb(t). (A.7) 

The physical meaning and conditions for the above hypothesis have been elaborated in Casimir 
(1945). 

Therefore, by (A. 2) we have 

E[a (g) (a(t + r) - a(t))] = E[ai ® ^[a(t + r) - a(t)|a(t)]] = r^[a ® Mh{t)] = rksM^ , 
E[{ai{t + r) - a(t)) (g) a] = E[E[ai{t + r) - a{t)\a{t)] a] = TE[Mh{t) a{t)] = rksM, 

and henceforth, by (A. 5), arrive at 

= M. (A.8) 



Applications to thermoelectric materials 

To apply the above general calculations to thermoelectric materials, we choose fi = (— V/i/T, —VT/T)'^ 
as the thermodynamic variables: a = fi, and assume that the system has unit volume. Taking the 
time derivative of (A.l) we obtain 

dS(a) 

where 7 is the entropy production rate per unit volume. From (2.21) we have 

jl = (jeJs) = -Sa, 
Therefore, if a and b are related by (A.6), then 

jl = -Sa = -SMb = -SMSa = -SMSfi. 
Comparing the above equation with (2.6) we have Li = —SMS, and henceforth 

Lf = Li 

since = M and S"^ = S. 
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